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Abstract 


The  frequency  distribution  of  a crystal  is  approximated  by 

combining  Van  Hove's  determination  of  its  analytical  nature  and 

2 

Mont  roll' a method  of  moments.  The  function  G(>i  ) is  represented 
by  an  expression  with  the  correct  behavior  at  the  singularities  and 
at  the  maximum  and  minimum  frequencies.  The  behavior  between  singular 
points  is  adjusted  smoothly  by  leaving  n undetermined  parameters. 

These  parameters  are  then  fixed  by  using  the  correct  first  r.  moments. 
As  a test,  this  procedure  was  applied  to  ''he  two  dimensional  square 
lattice  with  nearest  and  next  nearest  neighbor  interactions,  solved 
exactly  for  a particular  case  by  Montroll.  The  approximated 
distribution  function  had  the  right  form  at  the  end  points,  contained 
terms  of  the  appropriate  logarithmic  form  and  a jump  function  (with 
known  coefficients).  It  also  included  Legendre  polynomials  with 
unknown  coefficients,  which  were  determined  by  the  moments.  The 
difference  between  the  exact  and  approximate  distribution  functions 
was  a few  per  cent  using  only  the  zeroth  moment  (normalization). 

Using  higher  moments  produced  a gradual  increase  in  accuracy. 
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Int  reduction 


The  normal  modes  of  a lattice  are  eigen-vectors  of  a matrix  M 


ij 


denoted  simply  by  M,  where  i,j  are  indices  of  the  cell,  and  <*  , (£ 
are  indices  of  particles  in  a cell.  Using  the  translational  symmetry 
of  the  lattice,  the  problem  is  customarily  reduced  to  the  diagonalizaf 
of  a matrix  M (k),  denoted  briefly  by  M(k),  where  k is  the 

propagation  constant  of  the  wave,  and  the  index  (or  ^ ) takes  on 
^2  different  values,  where  0.  is  the  dimension  of  the  space  and  2 
the  n’unber  of  atoms  per  cell.  The  coefficients  of  M(k)  are  periodic 
functions  of  k,  with  the  periodicity  of  the  reciprocal  lattice. 

■v2 

The  normalized  density  of  eigenvalues  V associated  with  the 
matrix  M can  be  defined  by 


\ 


(1) 


G(  V ) = < Lr<  - M)/> 

where  the  average  of  a matrix,  ) , is  the  trace  of  that  matrix 

divided  by  its  dimensionality.  Using  the  uniform  spacing  of  the 
propagation  constant  over  one  Brillouin  zone,  eq.  (1)  can  be  reduced  to 


G(  S ) = 


where 


V -M(k))^  = trace 


f <(  ^ - M(k)>dk/  f dk  (2 

cf(V2-M(k )}/£?■  =(Jll) 


The  relation  between  the  frequency  distribution  function 


and 


the  distribution  function  of  the  squares  of  the  frequency  G(  V ),  is: 

g fv)  = XV  G(*x)  <:.) 

Since  the  integral  in  (2)  cannot  be  evaluated  explicitly  except 
12  3 1 

for  a few  special  cases  * ’ many  approximate  methods  have  been  use! 


p y.  Mcntrol  1 , .t.  Chem.  Phys.  it,  5 ('p  (19*+T). 

V.  A.  Bowers  and  H.  B.  Rosenstock,  J.  Chem.  Fhys.  l8,  1056  (1950). 

3 g.  A.  Newell,  J.  C'nem.  Phys.  21,  1877  (1953)- 

K.  L.  Rosenstock,  J.  Chem.  Phys.  21,  2061  (1953)- 
K R.  Rcsenstock  and  Gordon  F.  Newell  21,  l607  (1953)- 

, H.  £.  Rosenstock  and  H.  M.  Rcsenstock  21,  l6o8  (1953)* 

A di a cun si on  cf  some  of  the  methods  together  with  references  is  gi . 
by  A.  C.  Jlenzies,  Repts.  Prog.  Phys.  l6,  83  (1953). 
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to  study  Qt*)-  Blackman  approximates 

frequencies  at  a large  number  of  points 

finding  their  distribution.  Th*«3  required  a great  deal  of  labor,  whjon 

6 

has  to  be  repeated  for  every  set  cf  force  constants.  Houston  finds 
the  distribution  along  special  lines  in  reciprocal  space  and  interpolctf 
for  the  rest  of  the  Brillouin  zone.  This  method  introduces  some 
spurious  singularities  in  - We  shall  discuss  Montroll's  'moment 

method  in  some  detail  as  it  forms  part  of  the  basis  of  our  work. 

Using  (1)  the  nth  moment  of  G(x)  is  given  by 


^(yjby  calculating  the 
in  the  Brillouin  zone  and 


TT)n  = ^xnG(x)dx  =<Mn>  (4) 

or  taking  explicitly  the  average  over  the  propagation  constant  as  in  (2) 
'Mn  = ^ <Mn(k)>  dk/  (dk  (5) 

Montroll  was  the  first  to  recognize  that  the  moments , are  easily 
computed  from  traces  at  powers  of  M(k)  via  (5).  He  showed  that  the  first 
N moments  can  be  used  to  construct  the  best  approximate  distribution 
function  in  the  sense  of  least  squares,  if  the  approximating  function 
can  be  expressed  in  a power  series  or  a polynomial  series  contain  N 
terms.'*'1  For  these  two  choices  of  expansion,  the  least  squares  criterion 


•?  M.  Blackman,  Repts.  Prog.  Phys.  8,  11  (1941). 

6 W.  V.  Houston,  Rev.  Mod.  Phys.  20,  l6]  (1948). 

J T.  Nakamura,  Prog.  Theor.  Phys.  5,  213  (1950)* 

q E.  W.  Montroll,  J.  Chen.  Phys.  10,  21b  (1942). 

° E.  W.  Montroll,  J.  Chem.  Phys.  11,  48l  (1943)* 

4°  T.  H.  Walnut,  J.  Chem.  Phys.  22,  692  (1954),  calculates  the  raome-ts 
directly  from  M. 

4-  The  use  of  polynomials  as  expansion  functions  will  always  lead  t i 
the  same  result  as  the  use  of  a power  series,  for  the  same  number 
of  terns.  The  use  of  orthognal  polynomials  greatly  simplifies  tiie 
arithmetic. 
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is  equivalent  to  the  requirement  that  the  first  N moments  of  the 

12 

approximating  and  exact  distribution  agree.  The  rate  of  converge  oe 

of  such  an  expansion  is  determined  by  the  smoothness  of  the  exact 

function.  Hence  if  G(x)  is  smooth,  the  use  of  even  a few  moments  could 

be  expected  to  give  a good  approximation  .^However,  it  was  found  by 

Montroll^  when  he  did  an  exact  calculation  for  a two-dimensional 

square  lattice,  with  nearest  and  next  nearest  neighbor  interactions, 

that  the  frequency  distribution  function  was  not  smooth  but  contained 

13 

two  logarithmic  singularities.  Smollett  extended  Montroll's  result  to 

the  case  of  a two-dimensional  ionic  lattice,  taking  into  account  the 

long  range  Coulomb  forces  between  ions.  Simile”  results  were  obtained 

2 

by  Bowers  and  Rosenstock  for  the  frequency  distribution  of  vibrations 

lU 

perpendicular  to  the  plane  of  the  lattice.  Van  Hove  subsequently 

showed  that  these  singularities,  far  from  being  accidental,  are  a 

necessary  consequence  of  the  periodic  structure  of  the  lattice.  They 

occur  at  the  critical  points  cf^  , points  at  which  1 grad  J^k)l  = 0. 

In  two  dimensions  G(x)  will  have  at  least  one  logarithmic  singularity. 

In  three  dimensions  there  will  be  at  least  three  critical  points  x , 

c 

where  while  G(x)  remains  continuous,  its  derivative  G^(x  ) has  an 

c 

inverse  square  root  singularity. 


If  an  expansion  is  made  in  functions  oth^r  than  polynomials,  e.g.  a 
Fourier  series,  the  method  of  equating  moments  and  the  least  squ  res 
method  are  distinct  for  finite  N.  While  knowledge  of  the  first  n 
moments  determines  a Fourier  expansion  to  N terras,  if  the  mcmert 
conditions  are  applied,  this  knowledge  is  insufficient  to  lead  to 
a least  squares  solution. 


^ M.  Smollett,  J,  Chem.  Phys . 15,  575  (19^7). 
^ L.  Van  Hove,  Thys.  Rev  89,  11 89  (1953)- 


These  considerations  show  that  an  approximating  function  Ga(x)  which 

is  smooth,  such  as  a linear  combination  of  polynomials  which  has  the 

correct  first  N moments,  would  converge  slowly  to  C-(x)  with  increasi  .g 

N.  Hence  a large  number  of  moments,  and  a corresponding  large  nmo’  .t 

of  labor,  would  be  needed  to  get  a good  fit  to  G(x). 

We  decided,  therefore,  to  take  explicitly  into  account  the  analytic 

nature  of  G(x)  in  constructing  Ga(x).  Ga(x)  is  represented  by  an 

expression  with  the  correct  behavior  at  the  singularities  and  at  the 

end  points.  The  behavior  between  singular  points  is  adjusted  smoothly 

15 

by  leaving  n-undeterrained  parameters  to  be  adjusted  by  the  moments. 

In  this  way  the  moments  are  used  for  approximating  a smooth  function  so 

that  good  agreement  to  G(x)  might  be  expected,  even  when,  only  a small 

number  of  moments  are  used.  Essentially  the  same  procedure  was 

suggested  independently  by  Posenstock, ^ who  applied  it  to  the  body- 

centered  and  face-centered  cubic  lattices.  See  Section  III  for  further 

discussion  of  Rosenst^ck's  results. 

II.  Two  Dimensional  Square  Lattice 

As  a test,  this  procedure,  which  we  call  the  moment -singularity 

method  was  applied  to  a two-dimensional  . monatomic  square  lattice  with 

neighbor  and  next  neighbor  interactions.  The  distribution  function 

obtained  by  this  procedure  will  be  compared:  (a)  with  the  exact 

distribution  calculated  by  Mcntroll,^  (b)  with  the  distribution  obtained 

by  the  unmodified  method  of  moments. 

For  this  lattice  M(k)  is  a two  by  two  matrix1^ 

' (t-i)c1-tc1c2+i 

V tS1S2  l+(t-l)C2-tC]C2y 

-1 


M(k)  = 1/2 


ts1s2 


(O 


where  = cos  , S^  = sin  k^ , i = 1,2,  t = (l+«</2^  ) . and  )C 


IS 

The  fitting  between  singular  points  may  also  be  done  by  any  other 
method,  e.g.  Blackman's  which  yields  sufficiently  accurate  results 
from  the  singularities. 

H.  B.  Rosenstock,  On  Counting  Lattice  Frequencies,  Bull.  Amer.  PL 
Soc.,Vol.  29,  No.  4,  I 10.  2 

1 ' For  comparison  with  notation  in  reiercnce  l:ki  ->  3^,  x -»  f , 

G(x)  (2f)-l  (f)  we  measure  the  frequency  in  such  units  tnat 

the  largest  frequency  v{=  1. 

-U- 


me  neighbor  and  next  neighbor  forces  respectively.  The  two  r^cts  of 


the  secular  equation  are  (7) 

x±  = 1/4  ^2t(l-CpC2)+(l-t)(2-Cp-C2jt>[4t  S12S22+(l-t)2(C1-C2)2]  ^ 

The  nth  moment  of  G(x)  can  be  found  from  (5)  by  integrating  the  trace 

M“(k)  ever  the  Brillcuin  zone,  -n  kp,  kg  ^ n,  and  dividing  the  -result 
2 

by  2(2n)  . Thus,  the  first  moment  - , i6  (8) 

/ « r n 1 

71)  = (8n2 )■•*■]  dkp\  dk  2-(l-t)(cos  kp+cus  k2)-t  cos  kp*cus  k2i  =1/2 

1 7-tt  /-it  2 ( g " j 

Montroll  has  calculated  the  first  six  even  moments  /i2n  of  ^ ( V ) , in 

a slightly  different  way.  These  are  related  to  the  'Yfl  n by 

“2„  - L M2ng<  v =^oxnG(x)dx  =Wn  <>-) 

and  their  values  are 

1Y\0  = i,1Tlp  = 1/2,  TH2  = i/4  + (i-t)2  4-  t2  (xO) 

7f,3  = i/8  + 3 [ (l-t)2  + t2j  /16 

7T) u = 1/16  + (3/16)  f (l-t)2  +t2]  + (3/128)  f(l-t)2+t2]  + 

* — Y 

+ (5/128)  [ (l-t)2t2J 

'VTlc  = 1/1024(462  - 980  t+960  t2-265  t3+80  tU-5  t5) 

10 

An  explicit  form  for  all  higher  moments  was  found  by  T.  H.  walnut. 

Due  tc  the  symmetry  of  the  problem  it  is  sufficient  to  work  in 

one  quadrant  of  the  Brillouin  zone,  0 - k^,  kg  fr  n.  In  this  quadrant 

the  critical  points  of  x t (k)  are  at,  (kp,k^)  = (0,0),  (0,it),  (*,0), 

(it, it)  when  0 <- t 6 1/5.  (pp) 

The  values  of  x,  in  the  two  branches  x t (k),  at  these  critical  points 

are  (C,  1,  1,  1-tJ  and  (0,  t,  t,  l-t)  nespectively . For  t greater 

l8 

than  1/5  but  less  than  1/2,  the  x+  branch  has  one  additional  critic  1 
point  at  (cos  1(t-l/4t),  cos"1  (t-l/4t j]  and  the  value  cf  x+(k)  at  that 
point  is  (1  + 3t)2/(l6t). 

1®  We  shall  assume  as  did  Mcntroll  in  ref.  1 that  t £ l/2. 
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From  the  behavior  of  x+(k)  in  the  neighborhood  of  the  critical 
points  in  reciprocal  space  it  is  possible  to  deduce  the  form  of  G(x) 
near  the  critical  values  of  x,^9  including  the  value  of  G(0)  and  G(l). 

If  the  singular  part  of  G(x)  is  called  F(x),  then  it  can  be  shown  ttr  t 

(llo  . 

i B /.n  | x-t  | +Bo  #n  | x-(l-t)l  +B,H{x-(l-t)  ),  Otttl/5 
F(x)  = -i  1 3 

(C^nfx-tj  +C^  n I x-(  l-i-3t ) / l6t(  +C3H(x-(l-t)  ),  l/5^t-l,'2 

where  H is  the  Heaviside  unit  function  H(v)  = G,  y <a  0,  H(y)  = 1,  f . r 


y ;>  0 , and  (12b) 

Bj_  = -[t(l-2t)j  '^2/it2,  B2  = -(l-3t)  £t(l-t)(l-2t)(l-5t)  J’^Vn2 
B^  = -(2/(t«2)  )f  de-j  (l-3t)/t+(l+k2  si:i2d)1//2i  ’1,  k2=( 3t2-2t+l)/t2 
Cx  = - [t(l-2t)j'1/2/n2;  C2  = [l6t/«J(  1-t  j'  ["  3(5t-l)(3t+l)  j " ^ 2 

c3  = (i-3t)  [t(i-t)(i-2t)($t-i)J  _1/2 

The  value  of  G(x)  at  0 and  1 is 

G(0)  = (1+t)  Q t(l-2t-3t2)”]]  ’L’/2/n;  0(1)  = t‘1//2/jt  (12c) 

In  order  to  apply  the  moment-singularity  method  we  subtracted  F(x) 
from  G(x)  and  approximated  the  remainder  by  a linear  combination  of 
n Legendre  polynomials.  The  first  n-2  coefficients  of  the  expansion 
are  determined  by  the  requirement  that  first  n-2  moments  of  the 
approximate  function  agree  with  the  exact  cnes.  The  remaining  two 
coefficients  are  then  determined  by  fixing  the  end  points.  Hence, 

N 

when  N moments  are  known,  the  approximate  distribution  function  Ga  (::)  is: 


(N), 


N+l 


G '"'(x)  = F(x)  + T A.^P,  (2x-l),  0 £ x £ 1 
a kTO  K k 

where 

A^  - (2k+l)|  \ Pk(2x-l)G(x)dx-  ( Pk(2x-l)F(x)dx 

L '0  )0 


(-3a) 


(13b) 


For  a discussion  see  Section  III  of  this  paper,  also  Van  Hove, 
reference  (lU)  sections  II  and  III.  For  comparison  see  Montroll, 
reference  1. 
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for  k 4 N-l 


and 

A 


G(l)  - F(  l)  ••  f 


(S) 


. i . 

\ 


(11:) 


l 


-i  (-i)^VK'  - G(0)  - r(o)  - V1  (-1)1  A( 


(N) 


Note  tn 


it  Ay  ’ is  independent  ox  i'l  ch ? order  cf  t 


the  approximation. 


k(H)  _ ACC  ) 


= “k 


, for  k ^ N-l. 


when  the  unmodified  method  of  moments  is  used  to  approximate  G [ .-) , 

,(n) 


) , V 

the  approximation  function,  when  N noaonts  are  known,  G (x)  is 


.(N), 


^ in) 


. ’(x)  = > B F,  (rx-l) 


U0 


(m) 

the  By  are  again  given  By  eq.  (13b)  with  F(x)  set  equal  to  zero. 

A numerical  evaluation  of  the  and  the  By^  was  made  for  t = 1/3 

1 17 

This  is  the  va.lue  of  t for  which  Montroll  * evaluated  G(x)  in  a 
closed  fern,  so  that  comparisons  can  be  made.  This  will  indicate  how 
useful  the  moment -singularity  method  is. 


For  t = 1/3,  eqo.  (12)  and  (13)  yield 

GN(x)  = -3/n2  Sr.  I X-I/3I  -h/n2  /n  /x-3/4|  + £ A^N)P  (2x-l) 
a > 1 ' kT0  K k (1-) 

G(0)  = 2/rt ; G(l)  = O)1^/* 

(N.) 

The  values  of  the  Ay  in  different  orders  of  approximation  N * 6 are 

/ »T  \ 

summarized  in  table  1.  The  values  of  the  By  , which  are  independent  of 
N,  are  BQ  = 1,  B±  = 0,  B2  = -.U167,  B^  = 0,  B^  = -.2269,  B 5 = -.2290  (16) 
Comparisons  between  G(x),  G^(x),  G^(x)  are  presented  in  tables  2,  3,  and 

4. 

III.  Conclusicn 

As  can  be  seen  f rc. m Table  4 use  of  the  non-singular  method  of 
moments  results  in  a great  increase  in  the  accuracy  of  the  approximation 
ever  the  usual  method  of  nooervrs  for  the  same  number  of  moments. 
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The  additional  work  required  to  find  F(x)  was  small,  for  t = l/3 - 
For  the  particular  value  of  t,  x-(k)  can  be  expanded  in  a Taylor  se,  1 _> 
near  the  critical  points,  and  the  behavior  of  G(x)  near  the  critical 
points  can  be  read  off  directly  from  Van  Hove.  For  other  values  of  t, 
the  critical  points  (0,  0),  (jt,  n)  are  of  the  type  which  Van  Hove 
calls  generalized  critical  points.  The  behavior  of  x-(k)  near  these 
points  is 

x(k)  = xc  +|  f | Vp  ( 'f  / |||  ) + 0(  \tP  )»  % 

To  find  the  fora  of  G(x)  near  these  points  the  integral  in  (2)  has  to  be 

evaluated,  neglecting  terms  of  0(  Si?)-  Transforming  to  polar  coordinates 
in  reciprocal  space  the  integration  over  the  radial  variable,  can  be 
done  immediately  by  means  of  the  delta  function.  This  leaves  an 
integration  over  the  angle  variable  whose  evaluation  may  be  quite 

difficult,  if  it  cannot  be  f und  in  the  tables.  The  integral 

representing  in  eq.  (tit)  is  one  such  case. 

However  even  when  the  exact  fora  of  the  singular  part  of  G(x), 
such  as  the  coefficient  multiplying  the  logarithm  terms  in  two  dimensions, 
or  the  square  root  terms  in  three  dimensions,  is  not  known , converge-. ce 
would  be  improved  greatly  if  the  approximating  function  contained  terms 
of  the  correct  singular  form.  The  coefficients  of  these  terms  could 
then  be  determined  by  the  moments.  Thus,  for  the  case  considered  in 
this  paper,  it  would  have  been  preferable  if  the  approximating  function 
G (x)  consisted  of  polynomials  multiplying  the  logarithmic  terms  whose 
value  at  the  critical  points  is  the  correct  one.  This  w^uld  have 
taken  account  of  terms  in  the  distribution  function  of  the  form 

| x-x^t  . A smaller  number  of  moments  might  then  have  been  used 
to  obtain  the  same  accuracy.  This  was  not  done  because  of  the  extra 
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work  involved  in  solving  simultaneous  equations  for  the  coefficients 

of  the  multiplying  polynomials  since  no  orthognal  set  would  be  availu  le 

In  general  the  amount  of  work  required  to  find  more  moments  has  to  1 

balanced  against  the  work  involved  in  solving  simultaneous  equations. 

20 

There  is  a theorem  of  Feynman  which  is  useful  in  finding  some 
critical  points  without  diagonalizing  M(k).  This  theorem  states,  if 

A(°^)  V ( ) = \(  ^ ) V ( ^ ),  A(  )is  hermit  ian. 

Then 

ajuLil  = ( y aa.(3> 

^ A 2 v d ^ 

Thus  | grad  IS  (k)  I = 0 whenever  I grad  M(k)  I = 0.  This  can  happen 
only  at  those  k,  where  all  the  branches  have  critical  points 
s imu ltane  ou  s ly . 

A great  difficulty,  in  many  cases,  in>  finding  the  analytic  form 

of  the  distribution,  is  the  location  of  all  the  critical  points. 

Van  Hove's  arguments  predict  only  the  minimum  number  of  critical  points 

but  give  no  upperbound.  It  is  generally  easy  to  find  those  which  are  at 

symmetry  points  of  the  Brillouin  zone.  In  most  cases  investigated  thus 
1,2,3 

far  these  are  the  only  critical  points.  This  might  be  due  to  tue 

assumption  made  in  these  cases  that  the  forces  are  short  range.  The 
number  of  critical  points  might  be  expected  to  increase  with  the  range 
of  the  interaction  and  their  location  is  then  more  likely  to  be  at  non- 
symmetry points  of  the  zone.  Rosenstock^  has  devised  a method  for 

examining  the  presence  of  critical  points  inside  the  zone  from  the 
2 

behavior  of  V“(k)  on  the  boundary. 

For  three  dimensions,  with  short  range  forces,  when  the  critical 
points  are  at  the  symmetry  points  the  matrix  M(k)  is  easy  to  dlc-gonalize 
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at  those  points  to  find  the  critical  frequencies.  To  find  the  exact 
form  of  the  distribution  near  the  critical  points  degenerate 
perturbation  theory  has  to  be  used  to  find  the  eigenvalues  in  the 
vicinity  of  k,.  However,  since  in  three  dimensions  G(x)  does  not 
become  Infinite  at  xc  it  might  be  sufficient  to  use  the  right  form  of 
G(x)  near  x , the  coefficients  of  the  singular  part  would  then  be 
determined  by  the  moment  as  mentioned  above.  This  would  agree  with 
the  results  of  Rosenstock . ^ 
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Table  1 
N 

Value  of  the  coefficients  of  A£  of  the  Legendre  polynomials 
in  the  expansion  of  the  non-singular  part  of  distribution  function. 


»\k 

.N 

A0 

an 

A1 

N 

A2 

N 

a3 

N 

\ 

aN 

A5 

.N  N 

A6  A7 

GU) 

-.1306 

-.1599 

.1568 

g(2) 

- • 1306 

-.2620 

.1568 

.1021 

g(3) 

-.1306 

- . 2620 

.1409 

.1021 

.0159 

of1*’ 

-.1306 

-.2620 

. 1409 

.1388 

• 0159 

-.0368 

g<5) 

-.1306 

-.2620 

.1409 

.1388 

.0536 

-.0368 

-.0378 

g<6) 

-.1306 

-.2620 

.1409 

.1388 

.0536 

-.0429 

-.0378  .0061 

Table  2 


Convergence  of  the  moment-singularity  approximation  G^(  V^)  to  the 

& 

p 

exact  distribution  G(  \)  ) with  increasing  N.  (N  is  the  number  of 

moments  used.) 


2 

V 

| F( /) 
j 

GU)(V2) 

a 

g<2>(1,2) 

ci3)(  V2 

) G ( ) 

a * 

0<5>(  v2 

a 

(6)  2 2 

) G*  ( \J  ) G(  V ) 

O' 

• ^50f 

• 637 

.637 

.637 

.637 

.^37 

.637 

.C37 

.05 

.528 

.653 

.697 

.689 

.670 

.687 

.689 

.683 

• 1 

.620 

.689 

.763 

• 752 

.731 

.739 

• 739 

.760 

• 3 

I.36 

1.25 

1.34 

1.34 

1.36 

1.35 

1.35 

1.34 

1/3 

60 

CD 

CD 

CP 

CD 

ao> 

60 

00 

.4 

1.25 

1.08 

1.12 

1.14 

1.16 

1.17 

1.18 

1.17 

.5 

1.11 

• 897 

• 897 

• 911 

.911 

.937 

• 937 

.951 

.6 

1.17 

.939 

.890 

.901 

.879 

.891 

.887 

.908 

• 7 

1.52 

1.28 

i.20 

i .20 

1.17 

1.16 

1.16 

1.14 

• 75 

CP 

(P 

CP 

CD 

CD 

CO 

00 

CD 

• 9 

.942 

• 755 

.682 

.671 

.688 

.694 

.695 

• 733 

• 95 

.199 

.637 

• 599 

.585 

.604 

.621 

.619 

.63 

1 

.6853 

• 551 

• 551 

• 551 

• 551 

.551 

.551 

.551 

Table  3 


(N),  .2 

Convergence  of  the  unmodified  moment  approx imat ion  Gv<  ^ V ),  to  the 

2 

exact  distribution  G(  with  increasing  N.  (N  is  the  number  of 
moments  used.) 


2 

\)  j 

rv> 

)<£°<  v2>-43><  V2) 

)o<4>(  \Z>g<5)(v2)1 

^6)(  /))g<  >f) 

0 

i 

.583 

• 33d 

.103 

.637 

.05 

i 

• 7C2 

.655 

.646 

.683 

.1 

i 

.808 

.861 

• 770 

.740 

• 3 

i 

1.11 

1.13 

1.20 

1.34 

1/3 

1.14 

1.14 

1.21 

CO 

.4 

i 

1.18 

1.13 

1.20 

1.17 

• 5 

i 

1.21 

1.12 

1.12 

• 9rl 

.6 

i 

1.18 

1.13 

1.06 

.908 

• 7 

i 

1.11 

1.13 

1.07 

1.14 

• 75 

i 

1.05 

1.12 

1.09  ^ 

.9 

i 

.808 

.861 

.953 

. 73.: 

• 95 

i 

.702 

•655 

.664 

.631 

1 

i 

.583 

.332 

.561 

.551 

Table  4 
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